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Abstract 

At low energies the charge sector of one dimensional Mott insulators can be described in terms 
of a quantum Sine-Gordon model (SGM). Using exact results derived from integrability it is possible 
to determine dynamical properties like the frequency dependent optical conductivity. We compare 
the exact results to perturbation theory and renormalisation group calculations. We also discuss the 
application of our results to experiments on quasi-lD organic conductors. 
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1 ID Mott insulators 

The Mott metal-insulator transition is a paradigm for the importance of electron-electron interactions in 
condensed matter systems. It occurs in a variety of actual materials and has attracted much attention over 
the last fifty years |jl|. The underlying mechanism that drives the transition is by now well understood, 
but details on e.g. transport properties remain largely unknown in I? = 2 and D — 3 due to the lack 
of nonperturbative methods for treating strongly correlated electron systems. The situation is more 
fortunate in two cases: D — oo, where much progress has been made in recent years |^ and D = 1, where 
nonperturbative methods permit essentially a full solution of the problem. The ID case is the one we 
will be concerned with here. A full characterization of the Mott insulating phase requires the knowledge 
of dynamical correlation functions. The frequency dependent optical conductivity (t(w) is one of the 
most important examples from an experimental point of view. The behaviour of cr(a;) in the metallic 
regime is easily understood in terms of the Tomonaga-Luttinger theory Q . The situation in the Mott 
insulating phase is much more complicated due to the spectral gap that is dynamically generated by 
the electron-electron interactions. Here <j{uj) has until now only been studied by perturbative methods 
H, 1^, which break down in the most interesting regime of frequencies close to the optical gap. In these 
proceedings we use methods of integrable quantum field theory to determine cr(w) in ID Mott insulators 
for all frequencies much smaller than the bandwidth, which is the large scale in the field theory approach 
to the problem. Some of the results presented here have already appeared in . 
The paradigm of a ID Mott insulator is the Hubbard model 

H = -tY, {4,a^i+i,a + h.c.) +UY, nL]ni^i ■ (1) 

l-a I 

Here q^o- are fermionic annihilation operators of spin a at site I of a one-dimensional chain and 

"-i.cr — cl^ci,a- The ID Hubbard model is solvable by the Bethe Ansatz Q and the exact solution 
establishes the presence of a Mott transition at half-filling (one electron per site) and U = 0. In other 
words, for a half-filled band the model is insulating for any nonzero value of the on-site repulsion U, 
whereas it is metallic for a less than half-filled band. Many properties have been determined exactly 
1^, but the asymptotic behaviour of dynamical correlation functions are at present known only in the 
metallic |lOj and the gas phase jl^. In order to make progress in the insulating phase and to clearly 
expose the mechanism underlying the transition it is very useful to consider the scaling limit of the half- 
filled Hubbard model. This field theory limit corresponds to weak coupling and can be obtained directly 
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from the exact spectrum and scattering matrix of the lattice model [jlj, |l^ . It is defined by 

t^oo , U/t^O , M = — J — exp {-2TTt/U) fixed. (2) 

and has been studied by many authors p^ . On an operator level the field theory can be constructed 
along the lines of e.g. chapter 15 of One starts by splitting the electron operators into fast and slow 
components 

Q^o- — > y/oQ [ex.p(ikpx) Ra-{x) + e-}q){—ikpx) Lcr{x)] . (3) 

Here kp is the Fermi momentum (it is tt/2 for the half- filled band), and L„ are right and left moving 
electron fields and x = lao, where oq is the lattice spacing. Inserting this prescription into the Hamiltonian 
(|l|) one obtains 

n = j dx [L„ id^Ll - R„ id^Rl] +g J da; [l • I - J • j] 

+ 1 y [: I ■ I : + : I ■ I : - ; J ■ J ; - : J • J :] , (4) 

where vp = 2tao is the Fermi velocity and g — 2Uao- Here J and I arc the chiral components of SU(2) 
spin and pseudospin currents 

P = i^:i?ti?, :, = 



J' 
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^(i?|i?T-i?li?i) , J+=R\Ri ■ (5) 



Note that the Hamiltonian (|j) displays the required SO (4) symmetry [|l5[ of the half- filled Hubbard 
model. By employing the Sugawara construction, the Hamiltonian can now be split into two parts, 
corresponding to the spin and charge sectors respectively 

Ti. = Hc + Ti-s , 

He = [ dx[:I-I: + :i-i-]+g [ dxl-i, 



■Us 
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27rUs 



J dx[: 3 -J : + : 3 ■ J ■] - g J dx 3 -3 . (6) 



Here Vs = vp ~ Uao/2Tr and Vc — vp + C/ao/27r. Apart from the (marginally) irrelevant current-current 
interaction in the spin sector and the difference in spin and charge velocities, the Hamiltonian (^ is 
identical to the one of the SU(2) Thirring model (l6|. The latter is integrable Q and using its 
exact solution it is possible to determine dynamical correlation functions via the formfactor bootstrap 
approach. For example the optical conductivity was determined in this way and compared to numerical 
dynamical density matrix renormalisation group computations on the Hubbard model in [ p^ . Here we 
consider the optical conductivity for a more general model. We note that the optical conductivity is 
rather special in that for all cases considered here the electric current operator couples only to the charge 
sector, which greatly simplifies all calculations. Let us now consider an extended Hubbard model 



l;<7 
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where ni = ni^^ + ni^i- 



By repeating the above analysis we find that the scahng hmit takes the form 





= K + 
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j da; [: I • I : 4 


- : I • I :] 4 


- Jdx[gj^ {i+r + ri+) 




K 


2'KVs 
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j dx[.3-3: 


f : J • J :] 


— 2g± J dx 3 ■ 3 . 





where 



(8) 

g^ = {U-2V)a^, .g,, = (2[/ + 12y)ao , v'^ = v p + ^^^-^^^ . (9) 

Clearly spin-charge separation still holds, so that the two parts of the Hamiltonian (||) can be bosonized 
separately. Here we are interested in the charge sector only but note in passing that the spin sector is 
gapless if t/ > 2V and gapped otherwise. Applying the standard bosonization rules to in one 
arrives at the Hamiltonian of the SGM (|l^) ^ . 

The electric current operator of the lattice models dJ), (M) is given by 



-iet 
T 



[^■'^^J + l'^ " , (10) 



and does not commute with the above Hamiltonians. Using this becomes 

j^^Jdx[lHx)-P{x)] (11) 

in the field theory limit. From now on we drop the factor et/h which simply fixes the units in which we 
measure the current. 



2 The Sine-Gordon model 

The low energy physics of the charge sector of a general, pure, one-dimensional Mott insulator is described 
by the SGM as we have seen in the previous section for some specific examples. The action is given by 

SsG = J a; 1^(9,0)2 - 2/.cos(/3(/.)| . (12) 

Here we have chosen the normalization of the Bose field following |2^, |l] by specifying the short-distance 
behaviour of the two-point function as 

(^e^o.4>ix) g-»a,/,fo)^ ^ |a; - yp'^"' as \x - y\ ^ 0. (13) 
The SGM posesses a conserved (topological) charge 

where 

f - -^e^^^'d^d, (15) 

ZTT 

is the Noether current Q The electric current ( plj ) is proportional to the Noether current 



J = 5o0 , (16) 



^It is normalized such tliat soli tons and antisolitons have charges ±1. 
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where A is some nonuniversal constant ^. The SGM is integrable and has been studied in great detail 
over the last 25 years ^ Q . Let us review some results obtained from the exact solution (23| that 
we will need in the following. First of all some of the results obtained in the repulsive regime > 1/2 
appear to depend on the regularization scheme j24[ [ij] employed to deal with the UV divergences. Here 
we follow the method of |lq|, which is very natural from a field-theory point of view. 



The spectrum of the SGM depends on the value of the coupling constant /3 or, alternatively, on 



l-/?2 



(17) 



For < 1 the cosine term is relevant in the renormalization group (RG) sense and dynamically generates 
a spectral gap M in the excitation spectrum. In the repulsive regime 1 < ^ < oo the spectrum contains 
only charged particles of charge Q — ±1, which are called solitons and antisolitons. In this regime the 
spectral gap is related to the "optical gap" A i.e. the gap seen in the optical absorption spectra by 
A — 2M. At the so-called "Luther-Emery" (LE) point ^ = 1 the SGM is equivalent to a free massive 
Dirac fermion. In this limit the solitons become non-interacting particles and as we will see the Mott 
insulator turns into a conventional band insulator. In the limit ^ oo the sine-Gordon model acquires 
an SU(2) symmetry and describes the charge sector of the Hubbard model at half-filling in the limit of 
weak interactions as discussed above. In the attractive regime < ^ < 1 excitonic soliton-antisoliton 
bound states are formed and the spectrum becomes more complicated. Here we constrain ourselves to 
the repulsive regime and refer to [ p6[ for results on dynamical correlation functions in extended Hubbard 
models that correspond to the attractive regime in the SGM. 

As usual in a theory with relativistic dispersion e{p) — ^/p^ + it is useful to parametrize the 
spectrum in terms of a rapidity variable 9 defined by 



p = Msinh^, e — Af cosh^ 



(18) 



Let us distinguish solitons and antisolitons by an index e = ±. The exact 2-particle soliton-antisoliton 
scattering matrix is then given by p7| 



sinh ■ 



sinh ■^{9 — TTi) 

sinh ^ 
smh ^(9 — 7r^) 
a +b' ^a + b, 

roc sin(6lt/7r) sinh (^t) 



exp < —I 



t cosh (I) sinh (^^"j 



dt 



(19) 



The two-particle S-matrix (19) completely specifics all scattering processes in the SGM as multi-particle 
scattering is purely elastic and factorizes into two-particle processes. A convenient formalism for the de- 
scription of a dilute gas of particles with factorizable scattering is obtained in terms of the Zamolodchikov- 
Faddeev (ZF) algebra. The ZF algebra can be considered to be the logical extension of the algebra of 
creation and annihilation operators for free fermion or bosons to the case of interacting particles with 
factorizable scattering. The ZF algebra is usually introduced formally based on the knowledge of the 
exact spectrum and scattering matrix, which for the SGM was obtained in [^8[ For the SGM in the 
repulsive regime the ZF operators (and their hermitian conjugates) thus satisfy the following algebra 



h-92)Z'-{92)Z'^{9^) 



^We assume that A is nonuniversal because the electric current is not a conserved quantity for the lattice model. 
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Zl {e,)Zl {92) = z\ {B2)Z\ {Qx)siii [6^ - 62), 

Z-- i9,)Zl (92) = zl, (92)S^^ {92 - 9,)Z-'- {9,) + {2n)5ll5{9^ - 92), ei, £2 = (20) 

Here the two-particle scattering matrices Sl}'lf{9) are defined in Eg. ([l9|) and £j — ±. The factor 27r in 
the last equation stems from the normalization of the single particle asymptotic states (cf Eq.(p^)). 
Using the ZF generators a Fock space of states can be constructed as follows. The vacuum is defined 

by 

Ze{Om = . (21) 
Multiparticle states are then obtained by acting with strings of creation operators Zl{9) on the vacuum 

\9^...9,),_,,^Zl{9n)...Zl{eim. (22) 

We note that ( po| ) together with (^2|) implies that states with different orderings of two rapidities and 
indices Si are related in the following way 

\9n. . . 9k9k+l ■ ■ ■ 9i)e„...ekSk+i---ei = Se^^ettli^k - ^'fe+l)|6'ri • ■ • 0k+l9k ■ ■ ■ ^l) 6r,...e'^^-^e'^...ei ■ (23) 

The resolution of the identity is given by 

1 = E E r • • ■ (^1 . . . 0.\ . (24) 



3 Spectral representation of the optical conductivity 

An efficient method for the computation of correlation functions in integrable, massive quantum field 
theories is given by the form factor approach. This approach is based on the spectral representation, that 
expresses correlation functions in terms of an infinite series over multi-particles states. The two-point 
correlation function of some operator O can be written as 



{O{x,t)OH0,0))^Y.T. / (2\n ," exph^p,-x-e,-^ |(O|0(O,O)|g„...gi)e„...eJ^ (25) 
where pj and ej have the form ( |l8|) 



Pj = M smh9j, Cj — M cosh6'j, (26) 

and 

f{9^ . . . 0„)ei...e„ = (0|O(0, O)|0„ . . . (27) 

are the form factors (FF). They can be calculated using a set of "axioms" specifying their analytical 
properties (see ^ |3^ and H. Saleur's contribution to this volume). From a physical point of 

view we will be mainly interested in Fourier transforms of retarded two point functions. Their form factor 
expansions have the form 



OO /'OO 



X°{uJ,q)^i dx e*("+^'^)*-*«^([0(a;,i), 0^(0,0)]) 

J -OO Jo 

j 5{q~ MY. ^sinh 9 j) 5{q + MY^sinh 9 j) | 



— M J2j cosh 9j +ie ui + M J2j cosh 9j + ie 



(28) 
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We will be interested in the optical conductivity which is related to the imaginary part of the retarded 
current-current correlation function, ^■'(wjq), by 

a{uj > 0) ^ Im {xHlo, q ^ 0)} /lu. (29) 

Here the electric current operator j was defined in (|l^). Using (|2^) we obtain the following representation 
for the optical conductivity 

X S{mY^ sinh 9k)6{LJ - M ^ cosh 9k) 

k k 

= a2{uj) + (y^{uj) + ... (30) 

Here /■'(^i • ■ ■ ^n)ei...e„ = (0|j(0, 0)|6'„ . . . 6'i)£„...ej are the form factors of the electric current operator 
(p!^), (72 (w) and (74 (w) represent the contributions from 2 and 4-particle processes and the dots indicate 
processes involving higher numbers of (anti)solitons. We note that as a consequence of charge conjugation 
symmetry only intermediate states with an even number of particles contribute to this correlation function 
poj] . From ( po|) it is easy to see that only 2-particle processes contribute up to energies lj = 4Af, only 2 and 
4-particle processes up to a; = QM and so on. It has been previously observed for several models that the 
FF series converges much more rapidly than expected on the basis of such considerations ^ |36| . 

This may be understood in terms of phase space arguments [^3| , ^ . 

The n-particle form factors ( p7| ) for the current operator j'' in the SGM have been determined in ||3^ 
and can be used to calculate the first few terms in the expansion (|3^). The two particle FF is given by 



47r2Af\/]4,, cosh 



61+82 



/^(0i,02)+- = -neM-+ = — iT^id ; c(^i - (31) 

/? cosh ( («i^|)±l£j 



where 



C{0) =csinh6l/2 exp ( / dk -^^^ — I ^^^^ 

A; sinh ( sinh (fc) cosh (|) 

/4\3 fi /-oo sinh(|)sinh(i^A:) \ 1 

c=(l)exp U/ ; , ^ M ; rf^ -V- (33) 

\U \^4io fcsinh(f)cosh2(|)y SttCc 

The function C,{0) is analytic in the physical strip < Im^ < 2n and satisfies the following identities 

mSo{0) = a-O), C{0 - 2m) = C{~0) . (34) 
The four particle form-factor is far more complicated and can be represented as jsot] 

P{9,, 04)— ++ = -g-^MVAd^ H C(0fc - 0i) 

' k<l 

X J] (sinhp2+m-^„-i7r)/e])"' 

m,n— 1,2 

X 2 sinh[(6'4 + 6*3 - 6I1 - 6I2 - 27ri) /2(\ 

X exp(-i^e'fc) /" — ]J(^(a-6ifc)cosh(a- i^6lfc) 

X A(e2"/«|e2«i/«,e2»^/«|e2«^/«,e2«^/«)- (35) 
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Here C(^')j c, d have been previously defined, and 

(sin^ (M) sinh (^A:) + sinh^ 



= ^ 2 exp 



= ^ 2 exp 

2 f sinh' 
t— ^ 



i-oo 2 

' dfc- 



pOO 



dk- 



Jo 

inh^ (I) sinh (^fc) 
:sinh sinh (fc) 



fcsinh (^^^ sinh(fc) 



OO 

n 

n=0 



> 4 + 2 



n 



■|(n+l) 



(l + l(" + i)-^^)r(i + |n) 



(36) 



1 _ e'^^i/i 
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(37) 



The function (/3(^?) is even and, as can be easily see from the last expression, has poles and simple zeros 
in the physical strip, < ImO < 2tt, at the points 6 — i7r/2 + iTr^fc, k > and 9 — «37r/2 + iTr^fc, k > 1. 
As 9^ ±oo it behaves like ip{9) ~ 2exp(=F|(l + l/O^)- 

We note that the above integral expression is valid only in the regime ^ > 2. In general it is necessary 
to regularize the a-integral in ( ^sf ) but for ^ > 2 no such regularization is required. 

The form-factors with different ordcrings of £i, ...,£4 can be obtained from ( ^5| ) using the symmetry 
property (p3|). Using (^J) we obtain 



/^■(ei,e2,03,e4)-, + ,-,H 
6*3, ^2, ^4)-,-, 



^1, ^2, ^'a, ^4)-,-,+,+'S'^ 



- 1^2 



5+1(02 -^?3) 

(01, 02, 03, 04) + ,-,-,+ = 
P(02, 03, ^1, ^4)-,-, + ,+ - /■'(01, 02, 03, 04)-,-, + , + 5+l (01 - 03)^0(^1 - 02) 
5'l|(01- 02)5+1(01 -03) 
/■'•(01, 03, 02, 04)-,-, + .+ - P (01, 02, 03, 04)-,-, + , + 5^1 (02 - 03) g_ 
5Z+(0i- 02)5+1(02 - 03) 



(38) 



■+- (t/l - fc'2j- 



(39) 



The remaining orderings appearing in ( ^ ) can be obtained using the transformation properties of the 
current form factors under charge conjugation |30t] 



//i(01, 02n)ei,...,e2T. — ^/^i(^l, 02n)-ei, 



(40) 



We can now use the above expressions for the formfactors in the spectral representation for the optical 
conductivity (^o|). The two particle contribution is easily obtained by evaluating the 5- functions in (30) 



(72 (w) 



Here Q{x) is the Heaviside function, 

/(0) = /■'"(0/2, -0/2)+_ = -/^■(0/2, -0/2)_+ 

and 

= 2arccosh(a)) , w = lo/2M. 



(41) 

(42) 
(43) 
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As mentioned before, ([4l| ) is the full, exact expression for the optical conductivity for frequencies smaller 
than 4M. 

After some calculations the four particle contribution can be cast in the form 



where 



0-4 (w) 



e{Lu- AM) 
wl92 



cosh (?sinh 7 + lj^ — cosh 6* cosh 7 I — 1 



X [cosh^ 9 sinh^ 7 + w^l 



(44) 



arccosh((D — 1), b{6) — arccosh 



rj)2 - 1 - cosh^ e 



2 cosh 6* 



= In 



cosh(a/2) + exp(— 7) cosh^ 



2arccosh 



cosh^ 6 sinh^ 7 + a)^ — cosh 6 cosh 7 



The remaining integrals in (44) as well as the function ({9) have to be evaluated numerically. The 
latter is easily done to very high precision. The multiple integrals in the expression for the four-particle 
contribution are much more difficult to evaluate numerically. We estimate the precision of our results to 
be of the order of 10~^. 
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co/M 



Figure 1 : Two particle (solid line) versus one hundred times the four particle contribution (dashed line) 
in the form factor expansion for = 0.9. 

We have evaluated the two and four particle contributions to the optical conductivity for several 
values of The two and one hundred times the four-particle contributions for — 0.9 are presented 
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in Fig. 1 . Most importantly, the square root singularity, being a characteristic feature of band insulators, 
is suppressed by the momentum dependence of the soliton-antisoliton form factor and reappears only at 
the LE point = 1/2 (the behaviour of the optical conductivity in the vicinity of the LE point is shown 
in Fig. H and discussed in the next section). 

The rounding off of the singularity is due to the 9 dependence of the term (|4^). A similar behavior 
was previously noted for the Hubbard model at half-filling which corresponds to the special SU(2)- 
symmetric point /J^ = 1. We find that for anj^/3^ 7^ 1/2 there is a square root "shoulder" a{uj) (x \/uj — A 
for uj/A - K 1. 

The four particle contribution to a{Lo) is seen to be insignificant at low energies and becomes larger 
than the two particle contribution only at w « 180M for = 0.9. This suggests that the optical 
conductivity is well described by the combination of 2 and 4-particle contributions up to several hundred 
times the mass gap. Computation of higher order terms in ( ^0[ ) becomes cumbersome and probably of 
no physical interest, since the previous analysis suggests that they become important outside the region 
of applicability of the field theory approach to physical systems ^ Nevertheless it is interesting from a 
theoretical point of view to determine their importance at very high frequencies, which we will do using 
a different approach in section 0. 



4 Vicinity of the Luther-Emery point 

As we have indicated in the previous section, the LE point is quite special. Let us now discuss the 
behaviour of the optical conductivity in the vicinity of the LE point. This will exemplify some differences 
between Mott insulators and conventional band insulators. 

Let us recall that the SGM is equivalent to the Massive Thirring Model (MTM) |||| 

i^pj'^d^.ilj - I (-07^ ■0)^ - m-ipi^ (45) 



Smtm ^ I d X 
with the following identifications 
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TT 2 V ^ 



(46) 



Here tp — ^ ^^'^ ^ ~ tp^^^ are the usual two-component Fermi fields and the gamma matrices 
are chosen as 

j' = a\j'^2a', (47) 

where cr* are Pauli matrices. The fermions in the MTM correspond to solitons and antisolitons in the 
SGM. From ( p6| ) it follows that for ^ = 1 (/3^ = 1/2) the fermions become non interacting and the action 
( ]45| ) describes a two-band system of free, spinless fermions. In the absence of doping, in the ground 
state the upper band is empty while the lower one is completely filled. Thus the LE point describes 



a conventional band insulator. In fact, at ^ = 1 Eq. (41) yields a square root singularity above the 
threshold, in agreement with this interpretation. 

However, as soon as we deviate from the LE point this singularity immediately disappears and the 
optical conductivity vanishes at the optical gap. Close to the LE point (^) implies the following analytical 

•^Note that in order to determine behaviour just above the threshold we need to consider only the two-particle contribution 
to the optical conductivity. 

*The field theory approximation is appropriate as long as the frequencies considered are much less than the band width, 
which is the UV scale in the problem. In practical applications one is unlikely to encounter a situation where the band 
width is more than 1000 times the spectral gap. 
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expression valid for a) — 1 <g; 1 



cr{Lu) oc — 



- 1 



(48) 



The threshold behaviour of cr{uj) for several values of /3 is shown in Fig. 



100 




0.01 0.02 
(co-A)/A 



0.03 



Figure 2: Threshold behaviour of the optical conductivity close to the Luther-Emergy point for four 
different values of /3; /? = 0.72 (sohd), (3 = 0.73 (dotted), P = 0.74 (dashed) and /3 = 0.75 (long dashed) . 

We see that the square root singularity above uj — A for P'^ — 1/2 is replaced by a maximum occurring 
at ojmax = A + 0(72). As we take 0^ ^ 1/2 from above, uJmax approaches A and (j{ijJmax) diverges. 



5 Large energy behavior 

In order to obtain a rough estimate of the importance of the contributions involving six or more particles 
in the FF sum and to obtain a complete picture of the optical conductivity, we now calculate cr(w) at large 
energies. A convenient method to do this is Conformal Perturbation Theory (CPT)|^, Viewing the 
SGM as a Gaussian model 

Soauss = j CpX-^{d,4>f (49) 

perturbed by the relevant operator 2cos(/3(^), we can formally obtain correlation functions of the per- 
turbed theory by 

(e>(x.)o(x.,....,(x„)) . (50) 

/ —111. \ a-^x cos[b(p]\ 

(e '^J ^^')CFT 

We use a normalization in which 

nn I |4/5^| |4/3^ 

le ...e ...e )cft - n:U|xi-yjr ' ^ ^ 

The two-point function of electric currents can for example be obtained by considering 

(j(x) j(x')> = A{dr(k{x,T)dr'<l>{x',T')) = A lim ^-5,9,. (e'^^^^'^^e-'"^^^^')) . (52) 
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Using (pd) we have 



{j{^)]{^))^A\im—drdr'{e 



ia(p{x) — ia0(x') \ 



A 



I0(x') 



)CFT 



CFT 



(53) 



Eq.(|53|) provides a perturbative expansion in integer powers of the scale ^. Using the exact solution 
of the SGM it is possible to relate /i to the physical soliton mass M pOt] 



r(/32) 



7rr(l - /32) 



M 



V^r(i/2 + e/2) 



-| 2-2/3-^ 



2r(^/2) 



(54) 



In our case the CPT expansion is free of ultraviolet divergences (which would cause further compli- 
cations H, |l[ H) as long as f < oo but is known to suffer from infrared divergences. For example, the 
correlation function of bosonic exponents in the first line of ( |5^ ) develops infrared divergences at order 
0(^^") for 0^ < 1 — l/2n |4^. Here we only consider the term of second order in /i, which is free of 
divergences as long as 0^ > 1/2. For more general calculations one could follow the approach suggested 
in Q (see also [[43|). 

We can now calculate the current-current correlation function to leading order in CPT and use the 
result to determine the optical conductivity at large energies. The Wj-integrals in ( ^3|) are carried out 
using the methods of H|. We find 



Vr(2/32) 



t^r2(i-/32)r2(i + /32) 



r(l) 



-5) 



20Fr(i±£) M 



40^ 



(55) 



We emphasize that the ratio of the coefficients of the high- and low-energy asymptotics (^) , ( pl[ ) is Gxed 
[^,|21]. In other words, the amplitude of the power law in ( p5| ) is tied to the overall factor in ( ^ ) 
and the form factor expansion must approach the perturbative result in the large-w limit. A comparison 
between the form factor results and (^5|) is shown in Fig.^. We see that the asymptotic regime is not yet 
reached at energies as high as ^ lOOOM. In practical terms this implies that perturbation theory (PT) 
cannot be used to make contact with experiment. We note that the contributions due to intermediate 
states with 6,8... particles are all positive and will make the agreement of the form factor sum with PT 
in the region oj « lOOOAf only worse. 

A good way to overcome these deficiencies of bare PT is to carry out a renormalization-group (RG) 
improvement as performed in Following |^ we describe the SGM as the Sj^ perturbation of the 
SU{2)i Wess-Zumino-Wittcn model 



So 



'5(7(2) 



91 
27r 



^ / d^x Jo Jo 



J-J+ 



(56) 



Here Ja and Ja are the left and right currents of SU(2)i, normalized by the following operator product 
expansions 

1 

2^2 

1 



Mz)Jm = 7^ + o{i) 



Jo{z)J±{Q) = ±-J±{Q) + 0{l) 
J+{z)J_{z) = \ + -Jo(0) + 0(1). 



(57) 
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Figure 3: Comparison between the 2+4-particle contribution to the optical conductivity for 0^ — 0.9 
and the results from bare perturbation theory. 



The SU(2)-symmetric perturbation g±^ = g\\ corresponds to ^ ^ oo. The RG equations for the Sine- 
Gordon model can be cast in the form |20| 



dg^ _ g\\g± dg^ _ gj_ 
dt ~ 1 + ' dt ~ l + f 



where t is the RG scale. The solution of (58) is |20| 

1-/3' ^/9 



5J- = 4 



where 



/32 l-q 
(1_ 5)^2x2,3^-2 



> 511=2 



/32 1 _ 9 ' 



^ g(4-4/3^)(t-to)_ 



.4(1-/32), 

The perturbative result ( p5| ) is expressed in terms of f^y , g± noting that 

4 



(58) 



(59) 



(60) 



(61) 



and by considering (7 — > in (E^), (|60|). Finally one may fix the scale to by simply choosing it the same 
as in 1^ t — to — In "^t/ij^ j^ and obtain modulo terms of higher order in the couplings 



cr(w) 



2t^r2(2-/32)r2(i + /32) 



r(|)e^/^V? ' 
27/2r(i+£) 



4/3-^-4 



(62) 



^In order to choose to is a meaningful way one needs to calculate the first subleading term in the CPT expansion, which 
is outside the scope of these proceedings. 
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Figure 4: Comparison between the 2 and 2+4-particle contribution to the optical conductivity for 
/32 = 0.9 and the RG improved PT. 



The RG improved result ( p^ ) for a{ijj) is compared to the form factor result (or more precisely the 
sum of the two and four-particle contributions) in Fig.^ We see that on the level of accuracy of a log-log 
plot the agreement is rather good down to energies of the order of 5M. Combining the FF results with 
the RG improved perturbation theory we then obtain a good description of the optical conductivity over 
the whole frequency range. We believe that the good agreement of RG with the exact result even at 
rather small frequencies is probably a particular feature of the correlation function considered here. For 
other correlators like for example the spectral function there is no reason to believe that RG will work 
as well as it does for the optical conductivity. 



6 Applications 

Let us now turn to applications of our results for the optical conductivity. There are several materials 
which are believed to be one-dimcnsional Mott insulators in one of their phases. These are quasi-lD anti- 
ferromagnets like KCuF s [^ Carbon nanutubes possibly the striped phase in Lai.67Sro.33Ni04||47| 
and organic conductors |4§|, [49[ . Here we will concentrate on the latter due to the availability of extensive 
optical data. Of particular interest for our purposes are the (TMTTF)2X and (TMTSF)2X families, 
where X is an inorganic monoanion. These materials exhibit a rich phase diagram as a function of tem- 
perature and pressure and at sufficiently high temperatures are believed to be ID Mott insulators. 
The (TMTTF)2X family are presumably good examples of ID Mott insulators, but as the optical gap in 
these materials is of the order of the bandwidth a field theory description is inappropriate. 

For the (TMTSF)2X Bechgaard salts the ratio of optical gap to bandwidth is small and a field theory 
description is possible. The Bechgaard salts are highly anisotropic materials and can be modelled as 
weakly coupled, quarter-filled chains. It was suggested in [|j that at energies or temperatures sufficiently 
far above the 1D-3D crossover scale i?cr, the interchain coupling becomes ineffective and a description 
in terms of a purely ID model with charge sector (^2|) should be possible. This is a very nontrivial 
assertion as the microscopic lattice Hamiltonian appropriate for these systems is not given by a simple 
extended Hubbard model like (0), but includes an explicit dimerization [Q. The low-energy effective 
field theory obtained by bosonization is therefore not given by the simple form (|^), but contains other 
perturbing operators as well. At present there is also some uncertainty regarding the value of i?cr because 
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interactions can renormalize its bare value, set by the interchain couphng, downwards [pl| 




(0/M 



Figure 5: Comparison between the optical conductivity calculated in the SGM for 0^ = 0.9 (solid lines) 
and measured optical conductivity for (TMTSF)2PF6 from Ref. |^ (diamonds). The inset shows the 
same comparison on a logarithmic scale. 

There is a lot of ambiguity in fitting our results to the data. The value of the optical gap 2M is not 
known and, as discussed above, we cannot calculate the overall normalization of cr(w). We therefore use 
these as parameters in order to obtain a good fit at large lo (where the theory is expected to work best 
as 3D effects are unimportant) to the data for any given value of fi. We obtain reasonable agreement 
with the data for 0^ « 0.9, which corresponds to a Luttinger liquid parameter of Kp = /3^/4 « 0.23. This 
value is consistent with previous estimates (see the discussion in p8[). 

As is clear from Fig.^, the model (|l^) seems to apply well at high energies, but becomes inadequate 
at energies of the order of about 10 times the Mott gap (sa 1600/cm in (TMTSF)2PF6). Spectral weight 
is transferred to lower energies and physics beyond that of a pure ID Mott insulator emerges. 

There are at least two mechanisms that should be taken into account in this range of energies. Firstly, 
as mentioned before, a small dimerization occurs in the ID chains and will almost certainly affect the 
structure of (j(lo) around its maximum. Secondly, the interchain hopping is no longer negligible p2| and 
ought to be taken into account. 
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